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Abstract 

Previous  disagreements  concerning  a  linear  correlation  between  tbe  hydrogen 
(H)  concentration  and  the  extraordinary  refractive  index  Ug  in  proton-exchanged 
lithium  niobate  (LiNbOs)  have  been  resolved  by  partitioning  the  total  H  into  opti¬ 
cally  active  substitutional  H  and  optically  inactive  interstitial  H.  The  H  and  Li  spa¬ 
tial  variations  in  both  Z-cut  and  X-cut  crystals  were  determined  by  secondary-ion 
mass  spectrometry  (SIMS)  with  a  quantitative  evaluation  in  atoms/cm®.  These 
samples  were  proton  exchanged  in  neat  benzoic  acid  at  185°C  and  then  were  an¬ 
nealed  at  400“  C  for  times  t  from  6  to  240  min  in  wet  flowing  oxygen.  For  the 
Z-cnt  crystals,  fit  of  the  SIMS  measured  H  profiles  by  expressions  obtained  from 
the  diffusion  equation  for  diffusion  from  a  finite  layer  gave  a  substitutional  H 
diffusivity  of  Df  =  5.0  ±  0.3  x  10“^^  and  an  interstitial  H  diffusivity  of 

Df  =  1.4  ±  0.1  X  10“^^  cm^/s.  The  wet  flowing  oxygen  acts  as  a  constant  source 
of  interstitial  H  at  the  surface  with  the  diffusivity  Df  and  gives  an  integrated 
H  concentration  due  to  the  flowing  wet  oxygen  which  increases  as  \/t.  The  Li 
diffusivity  was  Df^  =  4.8  ±  0.2  x  10”^^  cm^/s  which  is  nearly  equal  to  Df.  For 
JA-cut  crystals,  the  substitutional  H  diffusivity  was  Df  =  3.4  ±0.2  x  10“^^  cm^/s 
and  the  interstitial  diffusivity  was  Df  =  1.3  ±  0.2  x  10~^^  cm^/s.  The  Ue  profiles 
were  evaluated  by  means  of  optical  prism-coupling  measurements  and  numerical 
simulations.  In  both  cases  of  crystal  orientation,  the  effective  index  diffusivity  is 
nearly  equal  to  the  diffusivity  of  substitutional  H.  Furthermore,  there  is  an  excel¬ 
lent  linear  relationship  between  the  profile  and  the  corresponding  substitutional 
H  distribution. 
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L  INTRODUCTION 

Crystalline  lithium  niobate  (LiNbOa)  is  an  important  material  for  modern 
optoelectronic  device  applications.  Due  to  its  high  electro-optic,  acousto-optic, 
and  non-linear  optical  coefficients,  this  material  has  found  widespread  use  in  op¬ 
tical  modulation,  optical  switching,  signal  processing,  and  frequency  conversion.^ 
Proton-exchange  (PE)  of  LiNbOa  is  an  effective  technique  for  formation  of  low 
loss,  optical  waveguides  in  this  material.  In  this  technique,  hydrogen  (H)  atoms 
in  the  liquid  melt  are  exchanged  with  Li  atoms  in  the  crystal  which  leads  to  a 
large  increase  in  the  extraordinary  refractive  index  However,  annealing  is 
typically  required  to  eliminate  compositional  instabifities  in  the  exchanged  layer 
and  to  restore  desirable  optical  properties  such  as  the  large  electro-optic  figure  of 
merit.^’^  In  addition,  the  depth  profile  of  Ue  in  the  PE-LiNbOa  layers  is  strongly 
dependent  upon  the  annealing  procedures.®  Consequently,  an  understanding  of 
the  relationship  between  the  redistribution  of  H  atoms  with  anneahng  and  the 
change  in  the  refractive  index  is  of  particular  importance  in  determining  optimum 
device  fabrication  procedures. 

Previously,  we  reported  on  the  correlation  of  H  depth  concentrations  with 
the  extraordinary  refractive  index  profiles  in  a  set  of  annealed  Z-cvA  PE-LiNbOa 
crystals.®*'^  The  H  distributions  were  determined  by  secondary  ion  mass  spectro¬ 
metry  (SIMS)  and  optical  characterization  of  the  refractive  index  was  done  through 
laser  prism-coupling  techniques.  It  was  found  that  the  measured  H  concentration 
decreased  more  slowly  with  depth  than  a  single  Gaussian  functioUj  which  is  the 
appropriate  solution  to  diffusion  from  a  thin  surface  layer.  However,  the  sum 
of  the  two  separate  Gaussians,  one  fit  to  the  low  H  concentration  region  and  an¬ 
other  to  the  high  H  concentration  region,  yielded  excellent  agreement  to  the  SIMS 
measured  profiles.®  In  analogy  with  diffusion  in  semiconductors,  the  larger  diffu- 
sivity  was  designated  the  interstitial  coefficient  Df  and  the  smaller  diffusivity. 
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the  substitutional  coefficient  Df .  It  was  possible  to  fit  the  effective  index  values 
for  the  observed  inodes  with  a  single  Gaussian  and  an  index  diffusivity  was 
determined.  The  average  value  of  was  equal  to  the  substitutional  H  diffusiv¬ 
ity  Z)f  and  the  index  depth  profiles  were  in  excellent  linear  agreement  with  the 
substitutional  H  distributions.^  These  results  suggest  that  the  change  in  the  ex¬ 
traordinary  refractive  index  depends  on  the  substitutional  H,  which  forms  crystal 
bonds,  rather  than  the  total  H  concentration  which  includes  interstitial  H. 

While  partitioning  of  a  diffusing  species  into  substitutional  and  interstitial 
components  is  common  in  semiconductor  materials,  this  concept  has  not  been 
adequately  addressed  in  the  annealing  of  PE-LiNbOa-  However,  there  is  experi¬ 
mental  evidence  that  two  H  components  are  present  in  PE-LiNbOa.  Incorporation 
of  H  on  substitutional  lattice  sites  as  well  as  on  interstitial  sites  has  been  reported 
by  Bolhnann  and  Stohr®  and  was  based  on  absorption  and  electrical  conduc¬ 
tivity  measurements.  Magnetic  resonance  measurements  by  Engelsberg  et  al.^° 
demonstrated  the  presence  of  both  substitutional  and  interstitial  H  in  annealed 
PE-LiNbOa  crystals.  Our  diffusion  analysis  indicates  a  similar  partitioning  of  H 
in  annealed  PE-LiNbOa  and  relates  the  substitutional  component  to  changes  in 
the  refractive  index. 

In  this  paper,  we  extend  our  prior  investigations  to  include  both  Z-cnt  and 
JC-cut  crystals  and  present  a  more  rigorous  model  for  analyzing  H  and  Li  diffu¬ 
sion  in  annealed  PE  LiNbOa  crystals.  For  this  model,  the  PE  region  represents  an 
extended  source  for  H  diffusion  into  a  semi-infinite  region  (substrate)  and  leads 
to  a  concentration  profile  which  is  the  sum  of  two  error  functions. The  Li  re¬ 
distribution  is  similarly  modeled  as  the  Li  diffusion  from  a  semi-infinite  source 
(substrate)  into  the  PE  region.  Again,  a  composite  fimction  which  is  the  sum  of 
two  error  functions  is  obtained.  While  a  single  diffusivity  leads  to  very  good  fits 
to  the  SIMS  measured  Li  profiles,  it  is  necessary  to  use  two  different  diffusivities 
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to  fit  tlie  H  profiles.  One  is  related  to  a  slower-diffusing  substitutional  component 
and  the  other  to  a  faster-diffusing  interstitial  component.  Furthermore,  we  find 
that  annealing  PE-LiNbOs  crystals  in  flowing  wet  oxygen  produces  an  additional 
H  source  which  needs  to  be  included  in  the  analysis.  The  refractive  index  profiles, 
that  are  determined  by  prism-coupling  experiments,  show  excellent  correlation 
with  the  substitutional  H  depth  profiles  in  both  Z-cut  and  X-cut  samples. 

Section  II  contains  details  of  the  experimental  procedures  used  in  this  investi¬ 
gation,  and  the  diffusion  models  for  analysis  of  H,  Li,  and  refractive  index  depth 
profiles  are  presented  in  Sec.  III.  In  Sec.  IV,  the  diffusion  analysis  is  apphed  to 
the  depth  profiles  found  in  annealed  Z-cut  PE-LiNbOa  samples.  The  diffusion 
analysis  for  the  depth  profiles  in  annealed  X-cut  PE-LiNbOa  samples  is  given  in 
Sec.  V.  A  summary  of  the  present  results  is  given  in  Sec.  VI,  including  a  discussion 
of  similarities  and  differences  found  in  the  analysis  of  Z-cut  and  X-cut  material. 
Questions  concerning  quantification  standards  in  the  SIMS  measurements  are  also 
addressed. 

II.  EXPERIMENTAL  PROCEDURES  AND  RESULTS 

Two  sets  of  LiNbOa  crystals  (Barr  and  Stroud)  were  used  in  this  study,  one 
Z-cut  and  the  other  X-cut.  A  proton-exchanged  (PE)  surface  layer  was  produced 
in  samples  of  each  set  by  inserting  the  samples  in  a  neat  benzoic  acid  melt  held 
at  a  temperature  of  185°C  for  70  min.  The  material  and  optical  properties  of  the 
PE- LiNbOa  samples  were  studied  in  their  as-exchanged  state  and  after  furnace 
anneahng.  The  samples  were  annealed  in  flowing  wet  oxygen  (-1.7 1/min)  at  400°C 
for  times  ranging  from  6  min  to  240  min.  Secondary  ion  mass  spectrometry  (SIMS) 
analysis  was  used  to  obtain  quantitative  H  and  Li  depth  distributions  for  these 
samples.  Refractive  index  profiles  in  the  exchanged  samples  were  determined  by 
optical  waveguiding  measurements. 
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A.  Hydrogen  and  lithium  depth  profiles 

Quantitative  atomic  depth  profiles  for  H  and  Li  in  the  as-exchanged  and  fur¬ 
nace  annealed  PE-LiNbOa  samples  were  annalyzed  using  a  Perkin-Elmer  6300 
quadrupole-based  instrument.  The  primary  Cs  ion  beam  energy  was  8  keV  and 
the  impact  angle  was  60°  measured  from  the  normal.  Lower  impact  energies  at 
higher  incidence  angles  have  been  shown  to  remove  ion  yield  artifacts  at  the  in¬ 
terface  between  the  PE  region  and  the  bulk  material.^^  With  this  method,  there 
are  no  “bumps”  in  either  the  H  or  Li  profiles  at  the  PE-exchange  interface.  No 
conductive  coating  was  applied  prior  to  analysis.  Profiles  of  H  were  acquired  by 
monitoring  negative  atomic  secondary  ions  whereas  the  Li  profiles  were  obtained 
simultaneously  by  monitoring  CsLi"^  molecular  ions.  Conversion  of  H  ion  count 
rates  to  concentrations  was  accomplished  by  using  relative  sensitivity  factors  cal¬ 
culated  from  analyses  of  a  H  ion-implanted  LiNbOs  standard  which  was  measured 
during  each  analyitcal  run.  The  H  profiles  were  background  subtracted  using  a 
background  determined  in  the  PE  samples.  The  Li  profiles  were  calibrated  by 
setting  the  value  of  the  profiles  in  the  unexchanged  region  equal  to  the  ideal 
atom  density  of  Li  in  bulk  LiNbOs  (1.89  x  10^^  atoms/cm~^).  Depth  scales  were 
calibrated  by  measuring  the  anal3rtical  crater  depths  with  a  stylus  profilometer. 

Figure  1  displays  the  H  and  Li  depth  profiles  that  were  obtained  for  the  Z-cnt 
PE-LiNbOs  crystals  studied  in  this  investigation.  These  samples  were  annealed 
in  flowing  wet  oxygen  at  400°  C  for  times  of  6,  15,  45,  90,  135,  and  180  min.  The 
H  depth  profiles  that  were  obtained  for  the  X-cut  samples  are  shown  in  Fig.  2. 

_  The  X- cut  samples  were  annea,led  under  the  same  conditions  for  45,  120,  and  240 
min. 

In  the  unannealed  PE-LiNbOg  samples,  the  H  depth  profiles  consist  of  a 
plateau  of  constant  H  atom  density,  followed  by  a  sharp  decrease  to  the  SIMS 
detection  limit  of  about  5  x  10^®  For  the  present  exchange  conditions, 

the  ^-cut  samples  have  an  exchange  depth  of  0.65  jum  and  a  peak  H  concentration 
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of  6  X  10^^  cm"®.  The  X-cut  samples  have  an  exchange  depth  of  0.75  jxm  and  a 
peak  H  concentration  of  3.5  x  10^^  cm“^.  While  a  step-like  H  profile  is  retained 
during  the  early  stages  of  anneahng,  particularly  at  relatively  low  temperatures 
of  approximately  275°  C,  annealing  at  higher  temperatures  of  near  400° C  rapidly 
produces  a  highly  graded  profile.®  The  graded  H  profile,  with  lower  concentration, 
extends  far  deeper  than  the  original  exchange  depth. 

The  behavior  of  the  Li  depth  distributions  is  somewhat  different  than  that  of 
the  H  profiles.  As  a  result  of  proton  exchange,  Li  is  depleted  from  the  surface 
region  of  the  LiNbOa  crystal.  The  Li  distribution  in  as-exchanged  crystals  is 
slightly  graded  at  the  surface  and  reaches  the  bulk  substrate  concentration  of 
=  1>89  X  10^^  cm“®  at  a  location  that  coincides  with  the  depth  of  the  PE 
plateau.^®  With  annealing,  Li  from  the  substrate  diffuses  into  the  PE  region  and 
the  Li  profiles  become  increasingly  shallow.  Eventually,  only  a  minor  difference 
exists  between  the  Li  density  near  the  surface  and  that  in  the  substrate  bulk. 

B.  Refractive  index  data 

Previous  investigations  have  indicated  that  the  extraordinary  refractive  index 
Tie  has  an  abrupt-step  profile  in  unannealed  PE-LiNbOs  samples.®’’^®  After  anneal¬ 
ing,  this  profile  becomes  highly  graded.  The  data  from  the  present  experiments 
support  such  behavior  for  rie.  The  PE-LiNbOa  crystals  were  examined  for  op¬ 
tical  waveguiding  using  a  prism-coupling  technique  with  a  HeNe  laser  operating 
at  A  =  0.6328  jwm.  Each  as-exchanged  sample  was  found  to  support  a  single 
TM  guided  mode,  which  is  consistent  with  the  exchange  depth  and  a  step-index 
change  of  approximately  0.10.  After  annealing,  the  samples  were  again 

characterized  using  a  prism  coupled  laser.  All  of  the  annealed  samples  functioned 
as  multimode  waveguides  and  up  to  five  modes  in  a  single  sample  were  observed. 
Effective  index  values  Ne  for  the  observed  modes  were  calculated  based  on  the 
prism- coupling  parameters.  An  inverse  WKB  method^^  was  used  to  determine 
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the  1/e  depths  d_i  of  the  modes  for  the  samples  that  supported  at  least  three 
propagatinjg  modes.  Depths  could  not  be  associated  with  effective  index  values 
for  as-exchanged  samples  and  for  samples  with  only  two  modes. 

The  optical  waveguiding  data  for  the  Z-cut  and  X-cut  samples  are  given  in 
Table  I  and  II,  respectively.  The  effective  index  values  Ng  for  the  observed  modes 
and  the  1/e  depths  d_i  which  were  determined  by  the  inverse  WKB  method  are 
listed  in  these  Tables. 

IIL  DIFFUSION  ANALYSIS  FOR  PE-LiNbOs 
A.  Hydrogen  diffusion 

The  diflPusion  of  H  from  the  PE  layer  is  governed  by  the  initial  conditions  of  a 
thin  surface  layer  with  a  constant  amount  of  diffusant,  as  shown  in  Fig.  3  (a),  and 
with  no  out  diffusion  of  H  from  the  crystal  during  annealing.  Mathematically,  the 
initial  condition  is  expressed  by: 

Cff(x,0)  =  C^  {otO<x<S,  (1) 


Off(x,0)  =  0 


for  X  >  6, 


where  is  the  initial  H  surface  concentration  after  proton  exchange.  With  the 
condition  of  zero  H  flow  from  the  surface,  the  boundary  condition  is  given  by: 

where  Cff{x,t)  is  the  H  concentration  in  cm“^  at  the  distance  x  and  time  i.  The 
resulting  diffusion  profiles  are  obtained  by  using  the  standard  diffusion  equation 
for  diffusion  from  a  finite  layer  into  a  semi-infinite  body  with  a  reflecting  boundary 
at  X  =  0  and  is  given  by:^® 


2VDi 


2VDi\' 
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where  erf  is  the  error  function,  x  is  the  distance  in  cm,  D  is  the  ditfusivity  in 
cm^/s,  t  is  the  diffusion  time  in  s,  6  is  the  exchange  depth  in  cm,  and  is  the 
H  surface  concentration  in  cm“^  at  t  =  0. 

When  S  <  y/Dt,  Eq.  (4)  reduces  to  the  Gaussian  function: 

qpe 

where  is  the  amount  of  H  per  cm^  in  the  PE  layer  and  is  given  by: 

Qh^  =  CffS.  (6) 

In  a  previous  paper,  a  Gaussian  function  was  used  to  analyze  the  H  depth  profiles 
in  ^-cut  annealed  PE  LiNbOs.®  When  a  single  Gaussian  function,  as  given  in  Eq. 
(5),  was  used  to  fit  the  H  profiles,  it  was  found  that  the  measured  H  concentration 
decreased  much  more  slowly  than  the  Gaussian  curve.  At  a  depth  of  approximately 
105,  the  calculated  Gaussian  curve  was  nearly  an  order  of  magnitude  less  than 
the  measured  profile. 

In  order  to  obtain  a  better  fit,  two  separate  Gaussian  functions  were  fit  to  the 
experimental  H  profiles.®  One  Gaussian  was  fit  to  the  low  H  concentration  tail 
region  and  a  second  Gaussian  was  fit  to  the  higher  H  concentration  surface  region 
so  that  the  sum  of  the  two  Gaussians  fit  the  measured  profile.  A  diffusivity  and 
surface  concentration  were  obtained  for  each  Gaussian  function.  The  diffusivity 
obtained  for  the  Gaussian,  which  was  fit  to  the  low  H  concentration  substrate 
region,  was  designated  the  interstitial  diffusivity  D,-,  and  the  diffusivity  for  the 
Gaussian  which  was  fit  to  the  high  H  concentration  surface  region  Was  designated 
the  substitutional  diffusivity  D,.  This  approach  led  to  very  good  fits  for  the 
experimental  H  data  and  indicated  that  two  species  of  H  are  involved  in  the 
diffusion  process. 

When  8  is  approximately  y/Dt,  it  is  more  appropriate  to  use  the  erf  represen¬ 
tation  given  in  Eq.  (4)  to  analyze  the  H  depth  profiles.  Figure  4  illustrates  the 
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difference  between  using  tbe  error  functions  of  Eq.  (4)  and  two  Gaussian  functions 
as  given  by  Eq.  (5)  to  model  tbe  H  depth  profile  for  the  Z-cnt  PE-LiNbOs  crystals 
annealed  for  6  and  45  min.  It  can  be  seen  that  while  the  overall  differences  are 
small,  they  are  more  significant  at  shorter  anneahng  times. 

As  with  the  Gaussian  function,  use  of  a  single  equation  of  the  form  of  Eq.  (4) 
did  not  produce  good  fits  to  the  experimental  H  data.  Vohra  et  al.^^  previously 
tried  using  such  an  expression  to  fit  H  diffusion  profiles.  Their  analysis  also  re¬ 
sulted  in  a  calculated  profile  which  decreased  less  rapidly  in  the  low  concentration 
tail  region  than  the  SIMS  measured  H  profile.  However,  very  good  fits  have  been 
obtained  to  the  data  shown  in  Figs.  1  and  2  when  two  such  expressions,  with 
different  diffusivities,  were  used. 

Since  the  function  which  is  given  in  Eq.  (4)  will  occur  often  in  this  analysis,  it 
is  convenient  to  represent  it  as  F(x,  S,  Dt): 

+  ^  (7) 

The  H  concentration  profile  GH{x^t)  is  then  given  by: 

CH(x,ty=C^;^,F(x,S,D,t)  +  C0F(xJ,Dit),  (8) 

where  and  are  the  surface  concentrations  of  the  H  diffusant  in  the 

substitutional  and  interstitial  components,  respectively.  The  sum  of  and 
is  the  total  H  concentration  introduced  into  the  crystal  due  to  PE.  The 

two  diffusivities  are  Ds  and  A* 

This  approach  implies  that  the  total  amount  of  H  in  the  annealed  PE-LiNbOs 
crystal  remains  constant.  However,  the  experimental  SIMS  data  indicated  that 
the  total  amount  of  H  actually  increased  with  annealing  time.  This  experimen¬ 
tal  result  was  foimd  with  both  Z-cut  and  X-cut  PE-LiNbOa  samples.  Since  all 
samples  were  annealed  in  flowing  wet  oxygen,  it  appears  that  the  flowing  gas 
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collected  water  molecules  and  introduced  H  into  the  samples  during  annealing.® 
Consequently,|Another  term  was  added  to  the  H  concentration  profile  Ch{x,  t)  to 
account  for  this  extra  source  of  H.  The  flowing  wet  oxygen  acts  as  a  constant 
source  of  H  at<the  surface: 

=  forOO,  (9) 

and  the  solution  of  the  diffusion  equation  from  a  constant  source  is  the  familiar:^® 

,  C0x,«)  =  Crerfc^,  (10) 

where  erfc  is  the  complementary  error  function  (1-prfV  In  gn.  (iQ),  is  the 

constant  H  surface  concentration  due  to  the  water  vaporJ%ls  assumed  that  the 

A 

H  due  to  wet  oxygen  annealing  enters  the  PE-LiNbOs  samples  primarily  as  an 
interstitial  species.  Therefore,  the  -complete  expression  for  the  H  concentration 
profile  is: 

Cjt{x.  t)  =  C'0F(x,  6,  Dst)  +  C^^F(x,  6,  Dit)  +  C^erfc^^.  (11) 

This  equation  is  the  analytical  expression  that  was  used  to  model  the  H  diffusion 
profiles  in  the  aimealed  PE-LiNbOs  samples.  Two  different  diffusion  mechanisms 
are  represent^:  a  finite  layer  diffusion  source  of  the  PE  region  and  a  constant 
surface  diffusion  source  from  the  water  vapor.  This  expression  also  involves  two 
species  of  H,  interstitial  H  and  substitutional  H,  each  having  a  different  diffusivity. 
The  interstitial  H  enters  the  samples  from  wet  oxygen  annealing  as  well  as  during 
PE.  The  substitutional  H  enters  the  crystal  only  through  the  PE. 

The  amount  of  H  diffusant  Qff(t)  per  cm^  in  the  crystal  after  annealing  for 
time  t  can  be  obtained  by  integrating  Eq.  (11)  from  z  =  0  to  x  -  oo.  There 
are  three  contributions  to  one  frona  the  substitutional  H  and  the  two  from 

the  interstitial  H.  Closed  form  expressions  can  be  obtained  for  the  integrals  of  the 
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collected  water  molecules  and  introduced  H  into  the  samples  during  annealing.® 
Consequently,  another  term  was  added  to  the  H  concentration  profile  C7^(a;,t)  to 
account  for  this  extra  source  of  H.  The  flowing  wet  oxygen  acts  as  a  constant 
source  of  H  at  the  surface: 

GH{^.t)  =  CW  fort>0,  (9) 

and  the  solution  of  the  diffusion  equation  from  a  constant  somce  is  the  familiar:^® 

CS:,(T,t}  =  (10) 

where  erfc  is  the  complementary  error  function  (l-erf).  In  Eq.  (10),  Cff"  is  the 
constant  H  surface  concentration  due  to  the  water  vapor.  It  is  assumed  that  the 
H  due  to  wet  oxygen  anneahng  enters  the  PE-LiNbOa  samples  primarily  as  an 
interstitial  species.  Therefore,  the  complete  expression  for  the  H  concentration 
profile  is: 

C„(x,t)  =  CS^F{x,S,D,t)  +  CS^F(x,S,J)<t)  +  Cff'edc^^.  (11) 

This  equation  is  the  analytical  expression  that  was  used  to  model  the  H  diffusion 
profiles  in  the  annealed  PE-LiNbOa  samples.  Two  different  diffusion  mechanisms 
are  represented:  a  finite  layer  diffusion  source  of  the  PE  region  and  a  constant 
surface  diffusion  source  from  the  water  vapor.  This  expression  also  involves  two 
species  of  H,  interstitial  H  and  substitutional  H,  each  haying  .a.differeiit  diffusivity, 
The  interstitial  H  enters  the  samples  from  wet  oxygen  annealing  as  well  as  during 
PE.  The  substitutional  H  enters  the  crystal  only  through  the  PE. 

The  amount  of  H  diffusant  Qnit)  per  cm^  in  the  crystal  after  anneahng  for 
time  t  can  be  obtained  by  integrating  Eq.  (11)  from  x  =  0  to  a:  =  cc.  There 
are  three  contributions  to  QH{i)i  one  from  the  substitutional  H  and  the  two  from 
the  interstitial  H.  Closed  form  expressions  can  be  obtained  for  the  integrals  of  the 
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error  functions: 

fOO 

J  F{xjS,Dt)dx  =  6,  (12) 

and 

Jo 

Equations  (12)  and  (13)  permit  writing  the  time  dependence  of  the  total  amount 
of  H  diffusant  as: 

QhW  =  QZ  +  OSf  +  CW^With.  (14) 

According  to  this  model,  the  total  H  diffusant  increases  with  annealing  time  as  the 
y/i  due  to  the  water  vapor  source.  The  H  contribution  from  the  proton-exchange 
remains  constant. 

B.  Lithium  diffusion 

In  order  to  analyze  the  Li  profiles,  the  Li  concentration  in  the  as-exchanged 
sample  was  assumed  to  be  a  step  function  as  represented  in  Fig.  3  (b),  and  that 
during  anneahng  there  is  no  outdiffusion  of  Li  from  the  PE-LiNbOs  crystals. 
Mathematically,  the  initial  conditions  are: 

Cx,(a:,  0)  =  (Tl,-  for  0  <  x  <  ^,  (15) 

and 

C'L,(x,0)  =  C'if^  forx>6,  (16) 

where  (7£,-  is  the  Li  surface  concentration  after  PE,  but  before  anneahng,  and 


.  is  the  Li  concentration  in  bulk  material.  The  .  boundary,  condition  at  the  surface 
is: 


dCiiix,  t) 


dx 


=  0, 


(17) 


a:=0 


where  Cu{x,t)  is  the  Li  concentration  profile  after  anneahng  for  time  t.  These 
conditions  represent  diffusion  from  a  semi-infinite  extended  source  in  the  substrate 
for  X  >  ^  into  a  surface  region  from  0  <  x  <  ^  with  a  reflecting  boundary  at 
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a:  =  0.  During  annealing,  there  is  a  concentration  gradient  which,  results  in  the  Li 
diffusion.  The  solution  to  the  diffusion  equation  for  these  conditions  is  the  sum 
of  two  complementary  error  functions.^^’^'^  Simulation  of  the  Li  diffusion  profiles 
with  the  use  of  MathematicdJ^  to  rapidly  compare  calculated  Li  profiles  with 
the  SIMS  measured  Li  profiles  showed  that  the  Li  surface  concentration  could  be 
taken  as  zero  rather  than  as  the  measured  value  of  =  0.5  X  10^^  cm“^.  This 
simpHfication  permits  representing  the  Li  diffusion  profile  as, 


Cu{x,t) 


erfc 


(^  +  ^)  ,  f  (^-®) 


(18) 


where  Du  is  the  Li  diffusivity.  With  Eq.  (18),  good  fits  to  the  SIMS  measured  Li 
profiles  for  the  annealed  PE-LiNbOs  samples  were  obtained.  A  single  value  for  the 
Li  diffusivity  Du  was  sufficient  for  fitting  the  experimental  curves  and  anneahng 
in  water  vapor  did  not  lead  to  an  additional  term. 

C.  Refractive  index  depth  profiles 

Several  studies  have  shown  that  the  extraordinary  refractive  index  depth  profile 
in  as-exchanged  PE-LiNbOs  is  approximately  a  step  function.®’^®’^®  With  anneal¬ 
ing,  the  index  depth  profile  becomes  highly  graded.  While  these  studies  are  based 
on  indirect  experimental  techniques,  it  appears  as  if  the  refractive  index  were  a 
physical  quantity  undergoing  a  diffusion  process  with  annealing.  Actually,  the 
change  in  the  index  profile  is  a  consequence  of  the  H  and  Li  diffusion.  Neverthe¬ 
less,  it  is  possible  to  model  the  diffusion  of  the  refractive  index  when  annealed 
with  the  following  initial  conditions;  ‘  . .  -  . . 


ne(a;,0)  =  nf^ 

for  0  <  a:  <  5, 

(19) 

ne(a;,0)  =  n^“^'' 

for  X  >  S. 

(20) 
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The  boundary  condition  at  the  surface  is: 


dne{x,t) 


dx 


=  0. 


(21) 


a:=0 


Here,  is  the  extraordinary  refractive  index  of  the  unchanged  bulk  crystal  and 
is  the  index  of  the  PE  layer.  The  resulting  solution  to  the  diffusion  equation 
for  these  conditions  is:^^’^® 


n, 


(22) 


where  Anf^  =  (nf^  —  is  the  initial  index  change  for  the  as-exchanged 

crystal,  and  Dn  is  the  index  diffusivity. 

As  with  the  Li  concentrations,  it  is  possible  to  model  the  depth  profile  of  the 
extraordinary  refractive  index  in  annealed  PE-LiNbOa  by  the  above  expression 
with  a  single  diffusivity  Dn-  No  additional  source  terms  were  required  due  to 
annealing  in  water  vapor.  The  total  index  diffusant  Qn  can  be  found  by  integrating 
[ne(a;,t)  —  from  a:  =  0  to  x  =  cx).  Integration  of  Eq.  (22)  gives: 


Qn  =  (n; 


-  (.r,PE 
e 


bulk 


)S. 


(23) 


This  quantity  depends  only  on  the  PE  conditions  and  should  remain  constant 
with  annealing. 

IV.  ANALYSIS  OF  DATA  FOR  Z-CUT  PE-LiNbOs 
A.  Hydrogen  depth  profiles 

-  -  The  H  depth  profiles  shown  in  Figs.  1  and  2  were  analyzed  using  the  composite 
expression  for  the  H  concentration  Cnix,  t)  given  in  Eq.  (11)  .  Since  the  following 
discussion  refers  to  ^-cut  material,  we  have  added  superscripts  to  the  various 
quantities: 


X 


2^fm' 


(24) 
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As  mentioned  previously,  this  expression  for  the  H  concentration  involves  two 
different  diffusivities  and  two  different  source  terms.  Two  different  diffusivities 
were  required  since  good  fits  to  the  H  depth  profiles  could  not  be  obtained  with  a 
single  diffusion  coefficient.  An  example  of  this  situation  is  shown  in  Fig.  5  where 
the  H  concentration  represented  by  the  single  Gaussian  of  Eq.  (5)  is  nearly  an 
order  of  magnitude  lower  than  the  measured  profile  at  a  depth  of  105^.  The  need 
for  two  different  H  source  terms  arises  from  the  experimental  SIMS  data  in  which 
the  total  H  in  the  crystal  increases  with  annealing  time  in  flowing  wet  oxygen. 
Consequently,  a  constant  surface  source  term  due  to  the  water  vapor  had  to  be 
used  in  addition  to  the  extended  H  source  from  the  PE  region.  In  unannealed 
.^-cut  PE-LiNbOa  samples,  the  PE  term  source  consisted  of  a  plateau  of  thickness 

=  0.65  )um  with  a  constant  H  concentration  of  6  x  10^^  cm“^.  These 
two  quantities  yield  a  total  amount  of  H  diffusant  =  3.9  x  10^^ 

cm~^. 

The  three  parts  of  Eq.  (24)  [PE  substitutional,  PE  interstitial,  and  wv  inter¬ 
stitial]  are  independent  and  combine  to  create  the  total  fit  to  the  SIMS  profiles. 
Figure  6  illustrates  this  procedure  for  an  annealing  time  of  135  min.  This  proce¬ 
dure  was  implemented  by  use  of  curve  fitting  routines  found  in  Mathematical^ 
Since  the  interstitial  H  will  diffuse  faster  than  the  substitutional  H,  the  interstitial 
mechanism  wiU  be  the  main  contributor  deep  into  the  crystal.  Near  the  surface, 
the  substitutional  mechanism  wiU  be  more  prevalent.  In  Fig.  6,  the  dotted  line 
represents  the  PE  interstitial  part,  the  lower -dashed  line  represents  the  PE  sub¬ 
stitutional  part,  and  the  dot-dash  line  represents  the  contribution  due  to  the  wet 
oxygen.  It  can  be  seen  that  at  greater  depths,  the  substitutional  part  from  the  PE 
becomes  negligible  and  the  interstitial  mechanism  (PE  plus  wv)  is  almost  the  only 
contributor  to  the  H  SIMS  profile.  The  solid  line  in  Fig.  6  represents  the  SIMS 
profile  and  the  upper  dashed  line  is  the  sum  of  the  three  terms  in  Eq.  (24).  Similar 
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numerical  fits  were  obtained  for  the  H  SIMS  profiles  for  the  other  annealing  times 
and  are  shown  in  Fig.  7. 

The  five  coefficients  available  Df,  and  Z)f )  to  fit 

the  data  for  an  individual  time  are  not  unique.  Trade-offs  can  be  made  between 
coefficients  and  still  give  an  excellent  fit.  Uniqueness  for  the  coefficients  is  de¬ 
termined  when  all  times  are  considered  using  a  trial  and  error  process  with  the 
following  constraints:  (1)  AU  of  the  coefficients  should  be  relatively  constant  over 
the  different  anneahng  times.  The  two  concentrations,  and  0^2^^ ■,  rep¬ 
resent  the  initial  surface  concentrations  of  the  PE  layer  and  represents  the 

constant  surface  concentration  due  to  wet  oxygen.  (2)  The  sum  of  ^rid 

must  equal  the  initial  total  PE  H  surface  concentration.  The  values  for 
these  quantities  in  Z-cui  LiNbOa  are  summarized  in  Table  III.  The  values  for  the 
diffusivity  of  Li  are  included  and  will  be  discussed  in  Part  B  of  this  Section. 

The  best-fit  values  that  were  obtained  for  the  diffusivities  [Df  ~  5.0  x  10”^^ 
cm^/s  and  Df  ~  1.4  x  10“^^  cm^/s)  are  nearly  constant  over  the  full  range  of 
annealing  times  from  6  to  180  min.  Since  Df  ~  ^Df ,  the  larger  diffusivity  was 
designated  the  interstitial  H  diffusivity  Df  and  the  smaller  diffusivity  was  desig¬ 
nated  the  substitutional  H  diffusivity  Df.  This  designation  is  reasonable  because 
in  most  materials,  the  interstitial  diffusion  coefficient  is  considerably  larger  than 
the  substitutional  diffusion  coefficient.^®  The  uncertainties  in  the  various  parame¬ 
ters  given  in  Table  III  for  H  are:  ±2.0  x  10^°  cm“^  for  ±2.0  x  10^°  cm~^  for 

±0.5  x  1020  cm-3  for  x  lO'^i  cm^/s  for  ,  and  ±0.3  xlO'^^ 

cm^/s  for  Df . 

The  amount  of  H  diffusant  Qfjit)  present  in  the  Z-cui  annealed  PE  samples 
after  annealing  for  time  t  is  obtained  by  integrating  Eq.  (24)  from  x  =  0  to  a:  =  oo. 
There  are  three  contributions  to  Qfjit),  one  from  the  substitutional  H  and  two 
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from  the  interstitial  H: 


QhW  =  Oh:™  +  '31/®  +  .  (25) 

Each  of  these  quantities  in  Eq.  (25)  is  plotted  versus  y/i  in  Fig.  8  and  is  labeled 
separately.  The  total  H  diffusant  increases  with  the  annealing  time  according  to 
the  SIMS  experimental  data.  The  H  diffusant  due  to  PE  is  nearly  constant  with 
anneahng  time,  with  about  75%  in  the  substitutional  component  and  25%  in  the 
interstitial  component.  The  contribution  to  Qff{t)  from  the  water  vapor  increases 
with  annealing  time  as  y/i  and  adds  to  the  interstitial  component.  As  seen  in  the 
figure,  the  amount  of  H  due  tb  the  annealing  in  a  flowing  wet  oxygen  ambient  can 
become  quite  large  and  even  equal  to  the  amount  from  proton  exchange.  However, 
according  to  the  model  presented  here,  the  substitutional  H  is  the  most  important 
H  component  because  it  forms  crystal  bonds  and  alters  the  extraordinary  refractive 
index. 

B.  Lithium  profiles 

The  Li  depth  profiles  that  were  shown  in  Fig.  1  for  the  annealed  Z-cut  PE- 
LiNbOs  samples  were  analyzed  according  to  the  model  presented  in  Sec.  III-B. 
The  Li  concentration  ,t)  represented  in  Eq.  (18)  was  used  to  fit  the  SIMS 
measured  Li  profiles  for  the  Z-ciit  PE-LiNbOa  samples  annealed  for  6,  15,  45,  and 
90  min.  The  results  are  shown  in  Fig.  9.  The  soHd  lines  are  the  SIMS  measured 
Li  profiles  at  the  designated  annealing  times  and  the  dashed  lines  indicate  the  ex¬ 
pression  in  Eq.  (18).  With  a  bulk  Li  concentration  of  1.9  X  10^^  cm“®,  the  values 
for  the  Li  diffusivity  Df,-  obtained  by  this  procedure  were  given  in  Table  III.  The 
average  Li  diffusivity  for  the  different  anneahng  times  is  =  4.8  ±  0.2  x  10“^^ 
cm^/s  which  is  nearly  equal  to  .  This  equality  of  diffusivities  suggests  that  the 
substitutional  H  diffuses  on  Li  sites  during  anneahng. 
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C.  Refractive  index 

As  described  in  Sec.  II-B,  an  inverse  WKB  method^^  was  used  to  determine 
the  1/e  depths  d_i  of  the  modes  for  the  samples  that  supported  at  least  three 
propagating  modes.  The  solid  circles  in  Fig.  10  correspond  to  the  Ne  and  d_i 
data  given  in  Table  I.  The  extraordinary  refractive  index  profile  given  in  Eq.  (22) 
may  be  re-written  with  Eq.  (7)  as: 

nf(x,t)  =  n“  +  Anf-™i?(x,<^Dft),  (26) 

where  is  the  extraordinary  refractive  index  =  2.204  at  A  =  0.6328  )«m) 
of  the  unexchanged  crystal  and  is  an  effective  diffusivity  of  the  index.  The 
sohd  lines  in  Fig.  10  are  the  best-fit  curves  to  the  experimental  data  by  Eq.  (26) 
and  the  values  used  for  Anf and  are  given  in  Table  IV  with  Df  for 
comparison. 

Because  ~  Df  and  because  the  same  analytical  expression  fits  both  the 
substitutional  H  profiles  and  the  extraordinary  refractive  index  profiles,  a  linear 
relationship  should  exist  between  Cff^g{x,t)  and  [nf(x,t)  —  The  scahng 

factor  relating  these  two  quantities  is  defined  by: 

nf(x,f)-nt“'*  =  a^C'|_,(a:,t).  (27) 

With  Df  ~  £>f,  this  relationship  becomes: 

.  ^  _  (28) 

The  values  of  determined  in  this  manner  are  given  in  Table  V.  A  comparison 
between  the  refractive  index  profiles  (solid  lines)  and  the  substitutional  H  profiles 
(dashed  lines)  is  shown  in  Fig.  10.  Only  for  <  =  15  min  can  a  difference  in  the 
two  curves  be  observed.  Otherwise,  the  correlation  between  the  refractive  index 
and  substitutional  H  profiles  is  excellent  for  times  from  45  to  180  min.  Based  on 
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this  analysis,  the  average  value  of  is  1.85  x  10“^^  (H  atoms/cm®)“^  when  the 
value  at  t  =  15  min  is  ignored. 

This  analysis  demonstrates  that  it  is  possible  to  model  the  depth  profile  of  the 
extraordinary  refractive  index  in  annealed  PE-LiNbOa  by  use  of  a  single  diffusivity. 
The  index  profiles  are  in  excellent  linear  agreement  with  the  substitutional  H 
profiles  in  Sec.  IV- A. 

V.  ANALYSIS  OF  DATA  FOR  X-CUT  PE  LiNbOs 
A.  Hydrogen  depth  profiles 

The  H  depth  profiles  for  the  X-cut  samples,  shown  in  Fig.  2,  were  analyzed  in 
the  same  manner  with  the  following  expression  for  the  H  concentration: 


C§{x,t) 


Dft)  +  C§f^Fix,  S^,  Dft)  +  C^-^exfc 


X 


iJWt 


(29) 

In  unannealed  X-cut  samples,  the  PE  extended  source  consisted  of  a  plateau  of 
thickness  6^  =  0.75  /itm  with  a  constant  H  concentration  of  =  3.5  x  10^^ 
cm“^.  These  quantities  give  the  total  amount  of  H  in  the  PE  layer  as  = 


=  2.6  X  10^^  cm-2. 

The  basic  fitting  procedure,  which  was  illustrated  in  Fig.  6  for  the  H  profile 
in  the  Z-cut  sample  annealed  for  135  min,  was  used  to  fit  the  X-cut  profiles. 
Numerical  fits  with  the  parameters  given  in  Table  VI  were  made  to  the  SIMS 
measured  H  profiles  for  the  annealing  times  of  45,  90  and  240  min,  and  the  results 
are  shown  in  Fig.  11.  The  SIMS  mea,sured  H  profiles  are  the  soHd  lines  and  the 
best-fit  curves  are  the  dotted  lines.  The  two  concentrations,  and 

represent  the  irdtial  substitutional  and  interstitial  H  surface  concentrations  in  the 
PE  layer.  Values  for  the  diffusivities  and  Df,  and  the  water  vapor  surface 
concentration  C§~'^  are  the  other  parameters  given  in  Table  VI. 

The  best-fit  values  that  were  obtained  from  the  H  diffusivities  for  the  X-cut 
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PE-LiNbOa  crystals  are:  ~  3.4  x  10“^^  cmVs  and  ~  1.3  x  10“^^  cm^/s. 

The  diffusivities  agree  well  over  the  full  range  of  annealing  times  from  45  to  240 
min.  While  the  substitutional  H  diffusivity  is  smaller  than  in  Z-cut  crystals,  the 
interstitial  H  dilfusivity  is  nearly  equal. 

The  amount  of  H  diffusant  Q^it)  in  the  X-cut  annealed  PE  samples  after 
anneahng  for  time  t  is  given  by: 

Q^{t)  =  Qi-/^  +  +  CS-'~yJiDXt/l,  .  (30) 

Each  of  the  quantities  in  Eq.  (30)  is  plotted  versus  V?  in  Fig.  12  and  labeled 
separately.  The  total  H  diffusant  increases  with  the  annealing  time  according  to 
the  SIMS  experimental  data.  As  with  Z-cut  material,  the  H  diffusant  due  to 
proton-exchange  is  nearly  constant  with  annealing  time,  with  about  75%  in  the 
substitutional  component  and  25%  in  the  interstitial  component.  Likewise,  the 
contribution  from  the  water  vapor  increases  with  annealing  time  as  the  y/i  and 
forms  part  of  the  interstitial  component.  The  amount  of  H  due  to  the  anneahng 
in  flowing  wet  oxygen  is  equal  to  the  amount  from  proton-exchange  at  about  90 
min. 

An  attempt  was  made  to  analyze  the  Li  depth  profiles  for  the  X-cut  samples. 
However,  due  to  the  shallow  Li  depletion  profiles  in  these  samples  and  the  long 
anneahng  times,  consistent  results  could  not  be  obtained  for  the  Li  diffusivity 
B.  Refractive  index 

~  Values  of  the  propagation  coefficients,  Ne  and  d_i ,  for  the  X-cut  samples  listed 
in  Table  II,  were  fitted  with  Eq.  (22)  which  may  be  written  as: 

nJix,t)  =  n^'‘^  +  Anf-‘‘‘^F{x,S^,D;t),  (31) 

where  is  the  effective  index  diffusivity  for  X-cut  samples.  Very  good  fits  to 
the  data  points  were  achieved  for  the  annealed  samples  as  shown  in  Fig.  13.  The 
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solid  circles  correspond  to  the  Ne  and  d_i  data  points  and  the  solid  curves  are  the 
best-fits  from  Eq.  (31).  The  values  used  for  Anf  and  are  given  in  Table 
VII  with  for  comparison.  As  with  the  Z-cut  samples,  the  average  value  of 
(4.6  X  10“^^  cm^/s)  from  this  analysis  is  nearly  identical  to  the  substitutional  H 
diffusivity  which  was  given  in  Table  VI. 

Since  ~  ,  the  refractive  index  change  can  be  linearly  related  to  the 

substitutional  H  distribution  as: 

=  (32) 

where  is  the  scaling  factor  for  the  X-cut  samples.  The  values  of  or^  determined 
in  this  manner  are  given  in  Table  VIII  and  a  comparison  of  these  profiles  is  shown 
in  Fig.  13.  The  substitutional  H  profiles  are  indicated  by  the  dashed  curves  and  the 
refractive  index  profiles  are  the  solid  curves.  Very  good  correlation  between  the 
refractive  index  changes  and  substitutional  H  profiles  exists  for  anneahng  times 
from  45  to  240  min.  The  average  value  of  based  on  this  analysis  is  3.9  x  10~^^ 
(H  atoms/cm^)"’^. 

The  analysis  for  X-cut  samples  provides  further  evidence  that  the  extrax)rdi- 
nary  refractive  index  profile  in  annealed  PE-LiNbOa  can  be  modelled  by  using  a 
single  diffusivity  and  that  these  profiles  are  in  excellent  linear  agreement  with  the 
substitutional  H  profiles. 

VI.  DISCUSSION  AND  SUMMARY 

Our  results  demonstrate  that  a  significant  amount  of  information  can  be  ob¬ 
tained  about  the  properties  of  annealed  PE- LiNbOa  through  a  careful  diffusion 
analysis  of  the  SIMS  measured  H  and  Li  depth  profiles.  Use  of  the  standard  dif¬ 
fusion  equation  with  the  appropriate  boundary  and  initial  conditions  leads  to  the 
identification  of  substitutional  and  interstitial  H,  and  values  for  the  diffusivities 
for  substitutional  H,  interstitial  H,  and  Li.  Analysis  of  the  total  H  concentration 
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provides  evidence  that  wet  flowing  oxygen  is  an  additional  source  of  interstitial  H, 
as  has  been  shown  in  other  ferroelectric  crystals.^^  This  procedure  has  also  made 
it  possible  to  establish  a  linear  correlation  between  the  refractive  index  change  and 
the  substitutional  H  concentration  in  the  crystal.  Previous  attempts  to  correlate 
H  distributions  with  the  refractive  index  profiles  have  resulted  in  non-Hnear^^  or 
piece- wise  linear^^  relationships.  However,  since  not  all  of  the  H  meaisured  in  the 
PE-LiNbOa  crystal  is  optically  active,  contradictory  results  have  been  reported. 
Furthermore,  annealing  the  crystals  in  flowing  wet  oxygen  contributes  to  the  H 
content  in  the  crystal.  The  additional  source  of  H  complicates  efforts  for  finding 
a  valid  correlation.  Our  model  is  based  on  a  conceptually  simple  approach  which 
provides  a  consistent  framework  that  establishes  a  linear  relationship  between  the 
refractive  index  profiles  and  the  substitutional  H  depth  distributions.  In  addition, 
the  area  beneath  the  refractive  index  profiles,  ais  represented  in  Eq.  (23),  remains 
nearly  constant  with  annealing.  This  result  agrees  with  data  on  similar  annealed 
PE-LiNbOs  planar  waveguide  layers, and  with  the  values  of  the  substitutional 
H  diffusant  as  plotted  in  Figs.  8  and  12. 

The  scaling  factor  a  that  relates  the  index  change  to  the  substitutional  H 
concentration  in  the  crystal  depends  upon  the  calibration  standard  used  in  the 
SIMS  analysis.  In  these  experiments,  a  hydrogren  (^H)  ion-implanted  LiNbOs 
sample  was  used  as  the  calibration  standard.  Several  attempts  were  made^«gisg 

.  PE-LiNbOs  samples  implanted  with  deuterium  (^H).  For  unexplained  reasons, 

spurious  profiles  of  were  obtained  and  the  results  could  not  be  used.  Moreover, 
there  was  a  change  in  the  H  surface  concentration  when  the  same  PE-LiNbOa 
sample  was  SIMS  analyzed  on  different  days.  rr,^y  4;^ 

ambiont  eondifciops  in-the  SIMS^e^pmcnt.  However,  the  depth  profile  curves  were 
of  identical  shape.  When  the  as-exchanged  PE-LiNbOa  samples,  and  several  of  the 
annealed  PE-LiNbOa  samples  were  SIMS  analyzed  on  the  same  day,  the  H  surface 
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concentrations  for  Z-cui  and  JC-cut  samples  were  found  to  be  equal:  = 

Ch~^.  This  result  means  that  while  the  diifusivities  determined  in  Sec.  IV  and  V 
are  not  sensitive  to  the  day-to-day  SIMS  calibration,  the  absolute  values  for  the  H 
surface  concentrations  do  vary.  Recent  data  obtained  from  H  forward  scattering 
analysis  of  these  samples  indicate  a  larger  H  concentration  in  the  crystals  than 
that  shown  in  Figs.  1  and  2?^  Further  experiments  are  needed  to  clarify  this  issue. 

Assuming  that  (7^“°  ~  values  of  a  for  Z-cut  and  X-cut  samples 

given  in  Secs.  IV  and  V  need  to  be  adjusted.  Using  the  value  for  cx^  (1.85  x  10“^^ 
cm®)  as  the  basis,  becomes  approximately  2.4  X  10~^®  cm®.  The  difference  in 
the  scaling  factor  for  Z-cut  and  X-cut  crystals  is  due  to  the  larger  index  change 
that  was  found  in  the  X-cut  PE-LiNbOa  samples.  It  should  be  noted  that  while 
there  is  some  uncertainty  in  the  overall  H  concentrations  and  in  the  values  of  a, 
the  basic  conclusions  of  this  work  are  unaffected.  There  are  two  species  of  H, 
substitutional  H  and  interstitial  H,  and  a  linear  correlation  exists  between  the 
refractive  index  change  and  the  substitutional  H  in  annealed  PE-LiNbOa. 
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TABLE  1.  Elfective  index  values  Ng  for  the  observed  modes  in  Z-cui  PE-LiNbOs 
substrates  after  annealing  at  400° C  for  time  t  (min).  The  1/e  depths  d_i  were 
determined  by  an  inverse  WKB  method  and  are  given  for  samples  that  supported 
at  least  three  modes. 


Mode 

<  =  15 

Ne 

min 

d-i 

<  =  45 

Ng 

min 

d-i 

<  =  90 

Ng 

min 

d-i 

<  =  135 

Ng 

min 

d-i 

<  =  180 
Ng 

min 

d_i 

TMo 

2.237 

1.0 

2.222 

1.4 

2.216 

1.9 

2.213 

2.3 

2.213 

2.8 

TMi 

2.214 

1.7 

2.212 

2.6 

2.209 

3.3 

2.209 

4.1 

2.209 

4.4 

TM2 

2.204 

3.1 

2.206 

4.0 

2.206 

5.2 

2.206 

5.8 

2.207 

6.2 

TMs 

2.203 

7.0 

2.204 

7.7 

2.205 

8.8 

2.205 

8.9 

TM4 

2.204 

12.7 
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TABLE  II.  Effective  index  values  Ne  for  the  observed  modes  in  X-cut  PE-LiNbOa 
substrates  after  annealing  at  400° C  for  time  t  (min).  The  1/e  depths  d_i  were 
determined  by  an  inverse  WKB  method  and  are  given  for  samples  that  supported 
at  least  three  modes. 


Mode 

t  =  45  min 
iVe  d_i 

t  =  120  noin 
N,  d_i 

t  =  240  min 

iVe  d_i 

TMo 

2.255 

0.91 

2.233 

1.50 

2.223 

2.17 

TMi 

2.232 

1.64 

2.224 

2.68 

2.218 

3.87 

TMa 

2.218 

2.58 

2.216 

3.68 

2.215 

5.23 

TMs 

2.210 

4.07 

2.212 

5.26 

2.212 

6.96 

TM4 

2.208 

8.29 

2.210 

7.87 

2.210 

9.72 
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TABLE  III.  Coefficients  for  the  fit  of  Eq.  (24)  to  the  SIMS  mecisured  profiles  for  H 
and  for  the  fit  of  Eq.  (18)  to  the  SIMS  measured  profiles  for  Li.  These  parameters 
are  for  Z-cut  PE-LiNbOa. 


Time 

. 

/^Z—wv 

Df 

Df 

Dl 

(mins)  (cm“^) 

(cm-3) 

(cm-3) 

(cia^/s) 

(cm^/ s) 

(cm^/s) 

6 

4.8 

X 

1021 

0.8 

X 

1021 

0.9 

X 

1021 

5.0 

X 

10- 

-12 

1.3 

X 

10- 

-11 

4.8 

X 

10- 

-12 

15 

4.5 

X 

1021 

1.1 

X 

1021 

1.0 

X 

1021 

5.0 

X 

10- 

-12 

1.5 

X 

10- 

-11 

5.0 

X 

10- 

-12 

45 

4.2 

X 

1021 

1.4 

X 

1021 

0.9 

X 

1021 

5.0 

X 

10- 

-12 

1.5 

X 

10- 

-11 

4.8 

X 

10- 

-12 

90 

4.2 

X 

1021 

1.4 

X 

1021 

0.8 

X 

1021 

5.0 

X 

10- 

-12 

1.4 

X 

10- 

-11 

4.8 

X 

10- 

-12 

135 

4.2 

X 

1021 

1.4 

X 

1021 

0.9 

X 

1021 

4.9 

X 

10- 

-12 

1.4 

X 

10- 

-11 

180 

4.2 

X 

1021 

1.4 

X 

1021 

0.8 

X 

1021 

5.0 

X 

10- 

-12 

1.4 

X 

10- 

-11 
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Table  IV.  Values  of  the  change  of  the  extrciordinary  refractive  index  Anf“^^  to¬ 
gether  with  the  diffusivity  of  the  index  and  the  diffusivity  of  the  substitutional 

H  I>f  for  ^-cut  LiNbOs. 


Time 

(min) 

D'i 

(cm^/s) 

Di 

(cm2/s) 

15 

0.100 

4.0  X  10-12 

5.0  X  10-12 

45 

0.084 

5.0  X  10-12 

5.0  X  10-12 

90 

0.072 

5.0  X  10-12 

5.0  X  10-12 

135 

0.072 

5.0  X  10-12 

4.9  X  10-12 

180 

0.084 

5.0  X  10-12 

5.0  X  10-12 
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TABLE  V.  Summaxy  of  the  extraordinary  refractive  index  change  Anf~^^  and 
the  substitutional  H  concentration  which  gives  the  scaling  factor  for 

the  Z-cnt  LiNbOs. 


Time 

(min) 

^n^-PE 

. 

(cm-3) 

(cm3) 

15 

0.100 

4.5  X  lO^i 

2.22  X  10-23 

45 

0.084 

4.2  X  1021 

2.00  X  10-23 

90 

0.072 

4.2  X  1021 

1.71  X  10-23 

135 

0.072 

4.2  X  1021 

1.71  X  10-23 

180 

0.084 

3.0  X  1021 

2.00  X  10-23 
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TABLE  VI.  Coefficients  for  the  fit  of  Eq.  (29)  to  the  SIMS  measured  profiles  for 
H.  These  parameters  are  for  X-cut  LiNbOs. 


Time 

(min) 

(cm-3) 

rX-PE 

(cin“3) 

/orA  ^wv 

(cm-3) 

Df 

(cm2/ s) 

Df 

(cm2/ s) 

45 

3.0  X  1021 

0.6  X  1021 

0.9  X  1021 

3.4  X  10-12 

1.1  X  10-11 

120 

2.8  X  1021 

0.8  X  1021 

0.9  X  1021 

3.4  X  10-12 

1.2  X  10-11 

240 

2.8  X  1021 

0.8  X  1021 

0.8  X  1021 

3.5  X  10-12 

1.5  X  10-11 
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Table  VII.  Values  of  the  change  of  the  extraordinary  refractive  index 
together  with  the  diffusivity  of  the  index  and  the  dilFusivity  of  the  substitu¬ 
tional  H  for  X-cut  LiNbOa. 


Time 

(min) 

(cm^/s) 

(cm2/ s) 

45 

0.122 

3.3  X  10-12 

3.4  X  10-12 

120 

0.110 

3.3  X  10-12 

3.4  X  10-12 

240 

0.100 

3.5  X  10-12 

3.5  X  10-12 
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TABLE  VIII.  Summary  of  the  extraordinary  refractive  index  change  and 

the  substitutional  H  concentration  C§~^^  which  gives  the  scahng  factor  a^. 


Time 

(min) 

-■  -^X-PE 

(cin-3) 

(cm3) 

45 

0.122 

3.0  X  1021 

4.07  X  10-23 

120 

0.110 

2.8  X  1021 

3.99  X  10-23 

240 

0.100 

2.8  X  1021 

3.57  X  10-23 
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Figure  Captions 

Fig.  1.  The  hydrogen  and  lithium  concentration  depth  profiles  in  2'-cut  PE- 
LiNbOs  samples,  as-exchanged,  and  after  annealing  at  400°  C  for  the  indicated 
times.  These  concentrations  were  determined  by  SIMS  measurements. 

Fig.  2.  The  hydrogen  depth  profiles  in  X-cut  PE-LiNbOs  samples,  as-exchanged 
and  after  annealing  at  400°  C  for  the  indicated  times.  These  concentrations  were 
determined  by  SIMS  analysis. 

Fig.  3.  The  idealized  depth  profiles  for  the  hydrogen  and  lithium  concentrations 
after  proton  exchange,  (a)  Hydrogen  depth  profile  with  the  surface  concentration 
C%  from  the  surface  to  the  exchange  depth  8.  (b)  Lithium  depth  profile  with  the 
surface  concentration  to  the  exchange  depth  8  and  then  the  constant  bulk 
value 

Fig.  4.  Comparison  of  the  fit  of  the  sum  of  two  Gaussians  as  given  by  Eq.  (5) 
(dashed  curves)  with  the  profiles  given  by  the  sum  of  two  error  functions  as  given 
by  Eq.  (4).  The  exchange  depth  8  was  taken  as  0.65  jxm  with  one  diffusivity  taken 
as  5.0  X  10“^^  cm^/s  and  the  other  taken  as  1.4  x  10~^^  cm^/s.  The  diflPusion  times 
were  taken  as  t  =  6  min  and  t  =  45  min. 

Fig.  5.  Attempted  fit  of  a  single  Gaussian  (dashed  curve),  to  the  measured  dif¬ 
fusion  profile  (soKd  curve),  for  a  .^-cut  PE-LiNbOa  sample  annealed  at  400° C  for 
45  min. 

Fig.  6.  Best  fit  of  the  composite  function,  as  represented  by  Eq.  (24),  to  the  mea¬ 
sured  H  dilfusion  profile  for  the  Z-cut  PE-LiNbOa  samples  annealed  at  400°  G  for- 
1 35  min.  The  various  quantities  plotted  in  this  figure  are  identified  in  the  inset 
to  the  figure. 

Fig.  7.  Best  fit  of  the  composite  function,  as  represented  by  Eq.  (24),  to  the  mea¬ 
sured  H  diffusion  profiles  for  Z-cut  PE-LiNbOa  samples  annealed  at  400°  C  for  the 
indicated  times.  The  SIMS  measured  H  profiles  are  given  by  the  soHd  lines  and 
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the  composite  functions  are  indicated  by  the  dotted  lines. 

Fig.  8.  The  amount  of  H  diffusant  in  .^-cut  samples  as  a  function  of  the  square 
root  of  the  anneahng  times.  The  upper  dashed  curve  is  the  least-squares  fit  to  the 
total  H  concentration  measured  by  SIMS  for  the  different  annealing  times.  The 
lower  dashed  curve  is  the  fit  to  the  integrated  substitutional  H  concentrations 
from  PE,  and  the  dotted  line  is  the  fit  to  the  integrated  interstitial  H  concen¬ 
trations  from  PE.  The  dot-dashed  line  is  the  fit  to  the  integrated  interstitial  H 
concentrations  from  the  water  vapor. 

Fig.  9.  Fit  of  the  complementary  error  functions,  as  represented  by  Eq.  (18),  to  the 
SIMS  measured  Li  diffusion  profiles  for  the  Z-cui  PE-LiNbOs  samples  annealed 
at  400°  C  for  the  indicated  times.  The  sohd  curves  represent  the  SIMS  measured 
Li  concentrations  and  the  dashed  curves  are  the  best-fit  by  Eq.  (18). 

Fig.  10.  Comparison  of  the  extraordinary  refractive  index  profiles  to  the  substi¬ 
tutional  H  concentrations  for  .^-cut  PE  LiNbOa  samples  annealed  at  400°  C  for 
times  of  (a)  15,  (b)  45,  (c)  90,  (d)  135,  and  (e)  180  min.  The  solid  curves  are 
the  Anf  depth  profiles  and  the  dashed  curves  are  the  substitutional  H  profiles. 
The  sohd  circles  represent  the  extraordinary  refractive  indices  calculated  from  the 
prism- coupling  measurements. 

Fig.  11.  Best  fit  of  the  composite  function,  as  represented  by  Eq.  (29),  to  the 
measured  H  diffusion  profiles  for  X-cut  PE-LiNbOs  samples  annealed  at  400°  C 
for  times  of  45,  120,  and  240  min.  The  SIMS  measured  H  profiles  axe  given  by 
the  solid  lines  and  the  composite  functions  are  indicated  by  the  dotted  lines. 

Fig.  12.  The  amount  of  H  diffusant  in  X-cut  samples  as  a  function  of  the 
square  root  of  the  anneahng  times.  The  upper  dashed  curve  is  the  least-squares  fit 
to  the  total  H  concentration  measured  by  SIMS  for  the  different  anneahng  times. 
The  lower  dashed  curve  is  the  fit  to  the  integrated  substitutional  H  concentration 
from  PE,  and  the  dotted  line  is  the  fit  to  the  integrated  interstitial  H  concen- 
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trations  from  PE.  The  dot-dashed  line  is  the  fit  to  the  integrated  interstitial  H 
concentrations  from  the  water  vapor. 

Fig.  13.  Comparison  of  the  extraordinary  refractive  index  depth  profiles  to  the 
substitutional  H  concentrations  for  X-cut  PE  LiNbOa  samples  annealed  at  400°  C 
for  times  of  (a)  45,  (b)  120,  and  (c)  240  min.  The  sohd  curves  are  the  An^  depth 
profiles  and  the  dashed  curves  are  the  substitutional  H  profiles.  The  solid  circles 
represent  the  extraordinary  refractive  indices  calculated  from  the  prism-coupling 
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